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Recent experiments on vesicles formed from 
block copolymers with liquid-crystalline side- 
chains reveal a rich variety of vesicle morpholo- 
gies. The additional internal order ("struc- 
ture") developed by these self-assembled block 
copolymer vesicles can lead to significantly de- 
formed vesicles as a result of the delicate inter- 
play between two-dimensional ordering and vesi- 
cle shape. The inevitable topological defects in 
structured vesicles of spherical topology also play 
an essential role in controlling the final vesicle 
morphology. Here we develop a minimal theoret- 
ical model for the morphology of the membrane 
structure with internal nematic/smectic order. 
Using both analytic and numerical approaches, 
we show that the possible low free energy mor- 
phologies include nano-size cylindrical micelles 
(nano-fibers), faceted tetrahedral vesicles, and el- 
lipsoidal vesicles, as well as cylindrical vesicles. 
The tetrahedral vesicle is a particularly fascinat- 
ing example of a faceted liquid-crystalline mem- 
brane. Faceted liquid vesicles may lead to the de- 
sign of supra-molecular structures with tetrahe- 
dral symmetry and new classes of nano-carriers. 

Amphiphilic block copolymers in water, like natural 
phospholipids, can self-assemble into various monolayer 
or bilaycr structures, such as micelles and vesicles 
In particular, rod-coil block copolymers, with a flexible 
hydrophilic chain and one or more rod-like hydropho- 
bic blocks, exhibit a rich morphology of structures, and 
therefore have significant potential to advance fundamen- 
tal science and drive technological innovations 0-[l2|. 
Among these rod-coil block copolymers, we are especially 
interested in liquid crystalline (LC) block copolymers in 
which the hydrophobic block is a nematic or smectic liq- 
uid crystal polymer (T3l - [20| . The self-assembly process of 
LC block copolymer vesicles is not completely controlled 
by the energetic and entropic interactions between differ- 
ent parts of the polymers. The in-plane LC order and the 
associated defect structure also play very important roles 
in determining the intermediate and final shape of vesi- 
cles. The tailor-design of both material properties and 
vesicle morphology by controlling the molecular struc- 
tures of the block polymers is state-of-the-art research 



in the fields of polymer science, materials science and 
chemical engineering. 

Some of the structures formed by these LC side-chain 
block copolymers in aqueous solution are rather counter- 
intuitive, such as faceted vesicles, nanotubes and com- 
pact vesicles with tiny inner space [HI, [l^ . In all these 
structures, the in-plane smectic order is clearly visible 
under Cryo-TEM. In this article we develop a theoretical 
explanation of the geometric structures of vesicles with 
in-plane nematic or smectic order. We present a simple 
model free energy as a functional of both the membrane 
geometry and the in-plane nematic order. Using both an- 
alytic and numerical methods, we then analyze the low 
free energy morphologies in various parameter regimes. 

Focusing on their overall shape we first look at the 
model free energy of a self-assembled monolayer as a 
functional of their shape and nematic order parame- 
ters [HI,!!!: 



(1) 



Here K is the Frank constant in the one-constant approx- 
imation, while D denotes the covariant derivative. H is 
the mean curvature and Hq is the spontaneous curva- 
ture, which is determined by the asymmetry in the sizes 
of the hydrophobic and the hydrophilic parts of the LC 
block copolymers. We shall choose the normal vector 
of the monolayer to point from the hydrophobic side to 
the hydrophilic side. Therefore H > means that the 
hydrophilic side is bending outwards. 

All three parameters K, k, Hq depend on the chemical 
structures of the block copolymers as well as their inter- 
action with the solvent in a complicated way. Further- 
more, strictly speaking, a nematic membrane is locally 
anisotropic. Therefore its Frank free energy is charac- 
terized by two constants: one for splay {Ki) and one 
for bend (A'a). Likewise, the bending energy as well 
as the spontaneous curvature should also be generically 
anisotropic, characterized by three bending constants 
and three spontaneous curvature components. Such a 
model is characterized by 8 independent parameters and 
is extremely complicated to analyze. For the sake of sim- 
plicity, we shall focus on the greatly simplified toy model 
Eq. ([1}, which captures the essential physics of nematic 
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vesicles which is the competition between the extrinsic 
bending energy and the two-dimensional Frank free en- 
ergy. 

A more important, conceptual issue is the following: In 
what sense can the vesicle morphology be understood in 
terms of minimization of elastic free energy Eq. ([1])? As 
is well known, the formation of vesicles is a complicated 
nonequilibrium process. Whether a certain property of a 
vesicle is distributed according to Gibbs-Boltzmann de- 
pends on the relevant experimental time scale, and on the 
time scale at which the given property equilibrates. At 
the stage of vesicle formation, individual molecules on the 
membrane can diffuse quite efficiently. Motion of liquid 
crystalline defects, however, requires coherent movement 
of all polymers on the vesicle, and is usually very slow. 
Hence we expect that the vesicle morphology achieves a 
local thermal equilibrium, where the shape and LC order 
minimizes the elastic free energy Eq. ([IJ (with appropri- 
ate parameters corresponding to the physical conditions 
under which the self-assembly takes place), subject to 
global constraints of given vesicle topology and LC de- 
fects distribution. We shall then enumerate all possible 
vesicle topologies and compare these free energy min- 
ima. It is interesting to note that in recent experiments 
by Jia et. al. [l^, multiple vesicle topologies were often 
observed using a given preparation method, suggesting 
that kinetics of self-assembly also played an important 
role in the selection of vesicle morphology. 

Smectic vesicles can also be viewed as nematic vesicles 
with bending constant much larger than splay constant. 
On a membrane with in-plane smectic order, therefore, 
the bending deformation of the nematic director field 
should vanish everywhere. Mathematically this is equiv- 
alent to h ■ Dfi = 0, that is, the nematic director lo- 
cally follows the geodesies. This is always possible, for 
an arbitrary but prescribed membrane shape, except at 
the core of nematic disclinations. For these configura- 
tions, the Frank free energy becomes independent of the 
bending constant. Hence Eq. ([T]) is also a toy model for 
membranes with in-plane smectic order, with the under- 
standing that K is the splay constant and the nematic 
director strictly follows the local geodesies. 

Minimization of the Frank free energy in Eq. ([!]) re- 
quires that the covariant derivatives of the nematic di- 
rector field vanish everywhere on the surface. As is well 
known in differential geometry, this is possible only if the 
Gaussian curvature vanishes everywhere, i.e. the surface 
is a developable surface. The family of developable sur- 
faces includes planes, cylinders, cones, and tangent de- 
velopable surfaces [28| . On the other hand, minimization 
of the bending energy in Eq. ^ leads to a constant mean 
curvature Ho . It is clear that the only geometry minimiz- 
ing both terms in the free energy in Eq. ([1]) is a cylindrical 
monolayer with a given radius 1/Hq. In the recent ex- 
ample of Jia et. al. |13j . for example, where only aqueous 
solvent is present at the final stage of assembly, mono- 
layer cylinders with very small radius (nanofibers) are 
observed. The inner space of the cylinders is completely 
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FIG. 1: Top: A tetrahedron compatible with +1/2 disclina- 
tion on each vertex can be constructed using a parallelogram, 
by folding along the dashed lines. A constant nematic direc- 
tor field in the unfold parallelogram is shown by the array of 
parallel straight lines. After folding up, the angles with same 
color circle around the same vertex. Bottom: the tetrahedron 
obtained via folding. There is exactly one +1/2 disclination 
on each vertex of the tetrahedron. 



filled by the hydrophobic parts of the polymers. In order 
to form monolayer cylinders with larger radius, the in- 
ner space has to be filled by solvent (or other polymers) 
that are friendly to LC blocks. If there is only aqueous 
solvent, and if l/i?o is not small, monolayer cylindrical 
structures with favorable spontaneous curvatures cannot 
pack space and therefore the system should form certain 
kinds of bilayer structures, where two monolayers with 
opposite orientation stack together. 

The free energy of a bilayer membrane can be obtained 
by adding up the free energies for two monolayers on both 
sides of the bilayer: 



(2) 



We shall focus on the morphology of bilayers in the re- 
mainder of this article. 

Without considering the boundary effects, a flat bilayer 
with uniform nematic order clearly minimizes both terms 
in Eq. ([2]) . The energy cost associated with the boundary, 
however, increases with the system size, and exceeds that 
associated with a closed vesicle with nonzero curvature, 
for sufficiently large systems [1^ . Close vesicles therefore 
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must form for sufficiently large bilayer membranes. 

The morphology of a bilayer is controlled by the com- 
petition between the extrinsic bending energy and the 
Frank free energy. We shall first limit the discussion 
to closed vesicles of spherical topology. Since the total 
Gaussian curvature is nonvanishing, the system is frus- 
trated and the Frank free energy competes with the bend- 
ing energy. First consider the limiting case K ^ n. The 
dominant contribution to the total energy is then the 
bending energy: minimizing this leads to a round spheri- 
cal shape. For a more realistic model where the bending 
energy is not isotropic, however, the shape will reflect the 
anisotropy of the bending moduli, leading to ellipsoidal 
shapes. The exact form of the shape as a function of the 
bending moduli is rather difhcult to calculate, however, 
and will not be treated in this article. Ellipsoidal vesi- 
cles are frequently observed in the experiments of Jia et. 
al. IJ], with the smectic layers all perpendicular to the 
long axis of the ellipsoid. It can be inferred from this 
observation that the bending rigidity is higher along the 
nematic director than perpendicular to the director. 

Let us now consider the opposite regime where K ^ k. 
In this case, the system should first minimize the Frank 
free energy, which leads to developable surfaces with van- 
ishing Gaussian curvature everywhere. This is clearly not 
possible due to Gauss' Theorem Egregium, which states 
that the total integrated Gaussian curvature of a sur- 
face with spherical topology is a topological invariant 
and equal to 47r. There are faceted polyhedral surfaces, 
however, for which the Gaussian curvature vanishes ev- 
erywhere but at a discrete number of (singular) vertices. 
These vertices are the ideal locations for orientational de- 
fects of the LC order (misery loves company) [23l- The 
total defect strength on a closed surface is also a topolog- 
ical invariant, according to the Gauss-Bonnet theorem. 
For nematic and smectic orders, this theorem dictates 
that on a sphere (or any other surface with the same 
topology), there are three possibilities for the structure 
of defects: 1) four disclinations each with strength +1/2; 
2) two defects each with strength -1-1; 3) one strength +1 
defect and two strength -1-1/2 defects. Now one needs at 
least four points to span a non-degenerate polyhedron; a 
tetrahedron in the minimal case. We conclude that in the 
limiting case K ^ k, the ground state morphology of a 
vesicle with spherical topology is a faceted tetrahedron, 
with a strength 1/2 disclination located at each of the 
four corners. This structure is indeed observed in recent 
experiments [l^, [2^ , as well as in our simulation, to be 
discussed in detail below. 

It is important to note that not all tetrahedra support 
a suitable nematic defect configuration. To ensure that 
the director field has vanishing covariant derivative ev- 
erywhere except at the four vertices, but including the 
six edges, the sum of the three angles surrounding every 
vertex of the tetrahedron must be 180°. This imposes 
three constraints on the geometry of the tetrahedron. 
[29| Since the set of all tetrahedral shapes (up to scal- 
ing the overall size) forms a five dimensional space, we 



see that the set of all fixed-size tetrahedra with vanish- 
ing covariant derivative everywhere except the vertices 
forms a two dimensional manifold. Fig. [T] illustrates how 
these tetrahedra, together with a nematic director field 
with vanishing covariant derivative, can be constructed 
by folding a parallelogram. These tetrahedra have the 
special property that all four triangular facets are identi- 
cal. All these structured smectic vesicles have vanishing 
Frank free energy. The degeneracy is lifted by different 
bending energies. The total free energy of the system is 
given by the sum of the bending energies localized on the 
six edges and the defect core energies localized at the four 
vertices. It is rather easy to see that for a given total sur- 
face area, the regular tetrahedron has a minimal value for 
the sum of all edge lengths. Thus the ground state mor- 
phology of a smectic vesicle with spherical topology is a 
regular tetrahedron when the bending rigidity is vanish- 
ingly small. The transition between different shapes is 
probably extremely slow, however, as it requires coher- 
ent motion of all four nematic disclinations together with 
the overall smectic layer texture. 

The edges and corners cannot be infinitely sharp in a 
realistic system. They are rounded by either the mem- 
brane thickness, the core size of a nematic defect or the 
small bending rigidity k. Likewise, the bending energy on 
the edges must be finite. In a realistic self-assembly pro- 
cess, the bending energy may also be partially relieved by 
preferential aggregation of large polymers on the outside 
and smaller polymers on the inside of the membrane near 
the ridges and corners. Faceted surface structures were 
studied previously in large viral capsids [1^, , which 
are formed by crystalline packing of proteins. There the 
faceting is energetically favorable because it reduces the 
in-plane strain energy of the crystalline order formed by 
the constituent proteins. What we have shown here is 
that a similar faceting can also be driven by the Frank 
free energy of LC order, despite their liquid nature. 

Another candidate for a low free energy morphology 
is a long cylinder of double layers (nanotube), for which 
the Frank free energy also vanishes. The bending energy 
is approximately given by 



□tube — kA/o'^, 



(3) 



where a is the radius of the cylinder. The total bending 
free energy is therefore linear in the membrane area. The 
faceted tetrahedron, on the other hand, has the total free 
energy 



tetrahedron 



4kL/6, 



(4) 



where L is the length of ridges and b is the radius of cur- 
vature of rounded out ridges. Since the area of a tetra- 
hedron grows quadratically in L, it follows that the total 
bending energy for a tetrahedron scales as the square root 
of the membrane area. Large faceted tetrahedral vesi- 
cles thus have lower free energy than nanotubes. Both 
morpholo gies , however, have been observed experimen- 
tally (l^. Il5l. [20| . Selection of vesicle morphology is 
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FIG. 2: Morphology of nematic vesicles at different bending rigidities. Left: The local Frank (a-c, g) and bending 
(d-f, h) energy contour plots; for a case of A'l = K3 = 1, a & d, k = 100.0; h &i e, n — 0.3; c & f , k = 0.05; for a case of 
Ki 7^ Ks, g &i h, K — 0.04. The total Frank free energies (in units of fcsT )are 62.83, 59.06, 57.1, and 65.27, respectively. 
The normalized total bending energies (in units of fcsT) are 12.55, 16.0, 18.58, and 50.34, respectively. Right; Calculated 
asphericities of nematic vesicles as a function of the inverse bending rigidity of 1/k. The inverse triangles are for a case of 
Ki — Ks and the triangles are for a case of Ki 7^ Kzi^K^/ K\ ~ 2.0). The empty circles represent locations corresponding to 
the morphologies of a (& d), b (& e), c (& f), and g (& h). 



also affected by kinetics of self-assembly, as we discussed 
above. 

In order to quantitatively investigate the ground state 
morphology of nematic vesicles, we develop a lattice 
nematic model on a deformable surface with spheri- 
cal topology and perform energy minimization by the 
method of simulated annealing Monte Carlo (MC). De- 
tails of the discrctizcd form of the free energy, whose con- 
tinuum limit is given by eq. ([2]), can be found in ref. [l^] 
and the Methods section. 

The simulation results for nematic vesicles at various 
bending rigidities show remarkable morphological tran- 
sitions, as displayed in Fig. [5] As the bending rigidity 
K decreases, the vesicles with an isotropic Frank elas- 
tic constant undergo substantial shape deformation: 1) 
the spherical morphology is found to be stable at large k 
(Fig. [5^, d); 2) ridges connecting four defects develop as k 
becomes smaller than 1 (Fig.l^b, c); 3) a tctrahedral vesi- 
cle forms at a vanishingly small k = 0.05 (Fig. [2}:, f). The 
faceting transition occurs near k ~ 1. The stable mor- 
phologies are determined by a delicate balance between 
the in-plane Frank energy and the bending energy as the 
surface deforms away from round. Indeed, as k decreases 
from 100.0 to 0.05, the Frank free energy falls from 62.83 
to 57.1 at the expense of bending energy which increases 
from 12.55 to 18.58. The Frank energy is localized near 



the four defects, which consequently induce deformation 
around the vertices. We note that while we continu- 
ously vary the bending rigidity k ranging from O.Ol to 
lOO.OfcfiT, the isotropic Frank elastic constant K is set 
to IfcfiT. Therefore, our simulations are entirely consis- 
tent with our prediction that spherical vesicles are stable 
in the regime oi K ^ k, whereas the faceted tetrahedral 
vesicles become stable in the other extreme K k. 

We also explore the effect of anisotropy in the Frank 
elastic constants by studying the regime in which splay 
dominates over bend. The smectic regime, as noted ear- 
lier, corresponds to the limit K3 ^ Ki. In the current 
simulation, the anisotropy is estimated to be K3/K1 ~ 
2.0 (see supplementary information). Although the shape 
transition trends are qualitatively similar for both the 
one-Frank constant and anisotropic cases, the anisotropy 
leads to a more dramatic shape transition, resulting in 
a considerably more faceted tetrahedral vesicle at a very 
small K ~ 0.04, as displayed in Figs. [2^ and h. In fact, 
the splay dominant nematic texture enhances the faceting 
more than the isotropic case does. This is clearly un- 
derstood by considering two membranes which possess a 
-|-1 disclination defect with pure splay and pure bending 
nematic textures, respectively. The pure splay always 
decreases the Frank energy by buckling out-of plane, be- 
cause it allows the defect to escape into the third dimen- 
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sion and thus better align the nematic directors. On the 
other hand, such out-of-plane deformation of the pure 
bending does not alter the Frank energy and therefore, 
the faceting of pure bending membranes is not favorable 
upon deformation. Note that we are restricting ourselves 
here to the case of isotropic bending rigidity. 

More prominent shape changes for vesicles with the 
anisotropic Frank elastic constants are clearly confirmed 
from a quantitative measurement of the asphericity (i.e., 
degree of deviation from the reference unit sphere geom- 
etry), which is defined as follows: 

(Ai?^) _ 1 (Ra — Ricf)^ 

"^'^^ ' 

where Ra is the radial distance of vertex a, i?rcf is the 
radius of the reference unit sphere, and N is the total 
number of vertices. The asphericities are averaged over 
10 simulation runs for each k and plotted in Fig. [2] as a 
function of 1/k. The plot exhibits relatively large devi- 
ations from its average values, especially at low bending 
rigidity. This is mainly attributable to the differences 
in the asphericity between the three possible ground 
state morphologies. Although faceted tetrahedral vesi- 
cles are expected to be the ground state for large sys- 
tem sizes, we have observed in our simulation three dif- 
ferent ground state morphologies, presumably due to 
its finite system size: i) an ellipsoidal vesicle with two 
closely bounded disclination pairs; ii) a flattened (square 
cushion-shape) vesicle with four -1-1/2 disclinations lo- 
cated approximately in one plane; and iii) a tetrahedral 
vesicle with four well separated -1-1/2 defects (see supple- 
mentary information). These three morphologies seem 
degenerate as the differences in their total free energies 
are within 0.5%. These vesicle shapes can be viewed as 
the precursors of the extreme morphologies at k — >■ 0, 
such as long fibrous cylinders, double layer sheets, and 
sharply faceted tetrahedrons, respectively. Finally, we 
briefly compare our simulation results with our recent 
experimental observations in Fig. [31 

In conclusion, we have studied the fascinating mor- 
phology of nematic/smectic vesicles, such as the faceted 
tetrahedron, nanofibcrs, and ellipsoids. Our theoretical 
and numerical studies provide the fundamental under- 
standing of formation of these novel structured vesicles 
and elucidate the shape-selective mechanisms. It could 
also pave the way for formulating guiding principles in 
designing nanocarriers with specific shapes, particularly 
utilizing the two-dimensional nematic order and the topo- 
logical defects, which are ubiquitous in closed vesicles. 

Methods 

To implement a deformable lattice model with spherical 
topology, we first introduce a reference sphere and tes- 
sellate it with a triangular mesh along with 12 requisite 
5-disclinations. Afterward a dual lattice of the triangu- 
lar mesh is constructed, each dual site being the center of 
mass of each plaquette formed by the original triangular 
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FIG. 3: Comparison between experimental obser- 
vations (a— c) and computer simulations (d— f). Left: 
Experimental results, a, a tetrahedron-shape smectic vesi- 
cle [231 ; b, a fat tetrahedron-shape smectic vesicle [l^; c, an 
elhpsoidal smectic vesicle [T^]- Right: Simulation results for 
a case of Ki/Ki ?a 2.0; d, k = 0.04; e, k = 0.1; f , k = 0.5. 
The contour plots show the distribution of the local Frank 
free energy. 

lattice. The details of the lattice geometry are illustrated 
in the supplementary information. Let rha to be the unit 
vector normal to the plaquette a. A director fia and a 
projection operator = fiafia are defined on each dual 
site with a constraint that it must be perpendicular to 
the plaquette normal: Na ■ riia = 0. 

Let dap be the bond length connecting two neighboring 
dual sites a and /3, and Sap be the spanned by the 
bond a/3. The discretized Frank free energy is then given 
by 

Fprank = K S^f^d^Ti [{Np - N^f] (6) 

<al3> 

The bond lengths dajs and the areas Sap are introduced 
to ensure that the lattice model is invariant under retrian- 
gularization, a necessary requirement for modeling a fluid 
membrane. In the continuous limit, the retriangulariza- 
tion invariance reduces to the usual reparameterization 
invariance of a fluid membrane. 

The discretized bending energy is given by 

-Fbending = K ^ S'^Tr K^, (7) 

a 

where is the extrinsic curvature tensor at site a, 
whilst Sa is the area of the plaquette a. The curvature 
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tensor Kq of each plaquette a can be calculated from the 
following three equations: 

ell - -e^ -K -e^ 

where is the vector pointing from vertex a to ver- 
tex (3, and e||,^ and e^^ are its components parallel and 
perpendicular to the plaquette normal rha- 

In the MC simulations, the deformable surface consists 
of 300 vertices, corresponding to 596 directors in all. The 
initial shape of the surface is a unit sphere and the initial 
director orientations are random. Each MC sweep con- 
sists of trial attempts to rotate each director and to move 
each vertex. The acceptance or rejection of a MC trial 
is determined by the standard Metropolis algorithm. All 
vertices are allowed to move along the radial direction 



with the angular positions of the vertices fixed. In order 
to preserve the total area upon surface deformation, any 
vertex moves making a total area change larger than 1% 
are rejected. Finally, once the surface is deformed by ver- 
tex moves, the orientations of directors are corrected by 
projecting them on the newly deformed plaquette before 
the new free energy is calculated. 
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